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Abstract 

cn , We exhibit the method for obtaining non perturbative quark and gluonic vacuum 
condensates from a model truncation of QCD. The truncation allows for a phe- 

^SJ ' nomenological description of the quark-quark interaction in a framework which 

^ . maintains all global symmetries of QCD and allows an -^ expansion. Within this 

^j. \ approach the functional integration over the gluon fields can be performed and 

CN| ' therefore any gluonic vacuum observable can be expressed in terms of a quark op- 

CS . erator and the gluon propagator. As a special case we calculate the mixed quark 

r~^ I gluon condensate gsigG^^o'^'^ q). We investigate how the value depends on the form 

0^ ' of the model quark-quark interaction. A comparison with the results of quenched 

r~| ' lattice QCD, the instanton liquid model and QCD sum rules is drawn. 

bj'' Key words: non perturbative methods in QCD, global color model, 

^ . Dyson-Schwinger equation, -^ expansion. 

PACS: 24.85.-Fp, 12.38.Lg, 12?38.-t, ll.15.Pg 



The existence of finite vacuum condensates such as the quark condensate 
(qq), the mixed quark gluon condensate g s{q G fiv^^^" q) i the gluon condensate 
o^siGfj.uG'^'') or the four quark condensate {qVqqVq) reflects in a direct way the 
non perturbative structure of the QCD vacuum. Their determination within 
a certain approach provides therefore important information about its ability 
to describe sufficiently the physics of strong interaction at low and interme- 
diate energies. One major field of application of these condensates are QCD 
sum rules [1-3], which are based on the operator product expansion (OPE). 
Hereby the vacuum condensates serve as the basic input parameters in which 
all the non perturbative effects are incorporated whereas all the short distance 
effects are treated within the standard perturbative diagram technique. The 
knowledge of the vacuum values at least of the lower dimensional quark and 
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gluonic operators is therefore essential for the apphcabihty of the QCD sum 
rule method. 

In this paper we want to focus on the mixed condensate gs{qG^j_pa^^'^q) of 
dimension five, which is among the least well determined ones. Previous studies 
of this condensate include QCD sum rules themselves where it was treated 
as fit parameter in the analysis of heavy-light quark system spectra [3,4], 
quenched lattice QCD [5] and, very recently, the instanton liquid model [6]. 

It is the aim of this letter to consider vacuum condensates in general and in 
particular the mixed condensate gs{qG^j^a^^q) in the framework of a trunca- 
tion of QCD which is based on an effective quark-quark interaction. For this 
let us consider the QCD partition function for massless quarks in Euclidean 
space 

Zqcd = J VqVqVA e-^QCDk-5,A] (^) 

with the Euclidean action 

SQCBh ^,A]=J dx {0 - tg4)q + \Gl^GlJ] (2) 

where 

Eq.(l) can be rewritten as 



A", 



Zqcd= DqDqe-'^rie'^^'Ssqi^-'il (4) 



where 



gWb1= /p^gJ^-4-M.-M.^JM-My ^ ^5) 



The functional VV[j] can be formally expanded in the current j" which leads 
to an expansion in terms of gluon n-point functions. 



w[j] = \ J dx,dx,j';\{x,)D;\i\{x,,x,)3;i{x,) + 
1 

3! 
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The first nontrivial contribution arises from tlie gluon 2-point function 

Df,{x, y) = D;1{x -y)= I VAA;{x)Al{y) e-^QCD[o,o,A] _ (7) 



It sliould be noted tliat in tliis definition of D^^, does not include quark loops. 

Our model truncation consists in truncating the series (6) after n = 2, i.e. 
we include only the gluon 2-point function. This defines an effective model 
based on the bilocal quark-quark interaction kernel D'^^^{x — y). This model 
truncation respects all global symmetries of QCD. It has been used extensively 
in the literature and is known as the global color model (GCM) [7-11]. The 
main features of QCD, which are lost by this truncation are local SU(3) gauge 
invariance and renormalizability. 

The partition function of this truncation is given by 



givhi^^giy) 



or 



Zoom = /p,P,-DAe-{^<^-'-^)-J'^^-*°-'4 , 

Eqs. (8) and (9) are connected by the functional integration rule 



(9) 



(10) 



where an obvious short hand notation for the integrations has been used. 

The gluon 2-point function D"^(x— y) is treated as the model input parameter, 
which, as we will discuss later, is chosen to reproduce certain aspects of low 
energy hadronic physics. For convenience we will use the Feynman like gauge 
D^y{x — y) = 5^y5°'^D[x — y) from now on. 

The next step is to apply the standard bosonization procedure , which consists 
in rewriting the partition function in terms of bilocal meson-like integration 
variables and expanding about the the classical vacuum, i.e. the saddle point 
of the action [7,8,12]. 

The resulting expression for the partition function in terms of the bilocal field 



integration is 2^gcm = / I^Be "^[^l where the action is given by 

B%x,y)B%y,x) 



S[B] = -TrLn 



S 



-1 



dxdy- 



2g,W{x-y) 



'111 



and the quark inverse Green's function, S* ^, is defined as 
S-\x,y)^^Jix-y)+A'B\x,y) . 



(12) 



The quantity A^ arises from Fierz reordering the current-current interaction 
in (8) 



''"''■^-{■^"V^-^4),, 



(13) 



It is the direct product of Dirac, fiavor SU{3) and color matrices 



® (73^-71'°-) 



(14) 



The saddle-point of the action is defined as |§|bo ^^^ i^ given by 



B'oix -y) = glD{x - y)ii K'S^ix - y) 



(15) 



where So stands for 5'[i3o]. These configurations provide self-energy dressing 
of the quarks with "rainbow" gluons through the definition S(p) = A.^Bq{p) = 
if \A{j?') — I] + B{p'^). The self energy functions A and B are determined by 
the "rainbow" Dyson-Schwinger equation 



A{p') - 1 



P 



8 f dS 



3 J (2vr; 
16 f d*q 



glD{p-q) 



A{q^)q ■ p 



4^s 



q^A^{q^)+B^{q^ 



(16) 



^^^'^-yIM-^^'^^'-'^ 



B{q' 



g2A2(g2) + 52(^2 



This dressing comprises the notion of "constituent" quarks by providing a mass 
M{p^) = B{p'^)/A{p'^), refiecting a vacuum configuration with dynamically 
broken chiral symmetry. 



We will calculate the vacuum condensates from the above saddle-point expan- 
sion, that is, we will work at the mean field level. This is consistent with the 
lowest order -^ in the quark fields for a given model gluon 2-point function. 
It should be noted, however, that the gluon 2-point function D itself contains 
all powers of ■^. 

The mesonic modes can be obtained within in the standard RPA technique 
by considering small fiuctuations around the saddle point configuration and 
solving for the eigenmodes of the arising quadratic kernels. This leads then 
automatically to the ladder Bethe Salpeter equations for the qq bound states. 
Mesonic properties have been extensively studied in refs.[7,8]. In ref.[12] a de- 
tailed investigation of the low energy sector was performed by deriving the 
general form of the effective chiral action for the SU(3) Goldstone bosons 
and determining /tt and most of the chiral low energy coefficients (Gasser- 
Leutwyler coefficients) Lj, which, in turn, determines the physics of the tt, 
K and r] mesons at low energies [13]. This was done for various forms of 
the model gluon 2-point function D{q^) all of them reproducing f.,^. In the 
present approach we are following this spirit by choosing gluon 2-point func- 
tions which have this feature and moreover reproduce values for the chiral 
coefficient Lj, i = 1,2, 3, 5, 8 which are compatible with the phenomeno logi- 
cal ones [13]. 

It is now rather straightforward to calculate the vacuum expectation value of 
any quark operator of the form 

Qn = {qh^f^h(l^l) (^iaAj^U^a) • • • (gj„AJj„gi„) (17) 



in the mean field vacuum. Here the A*^*^ stands for an operator in Dirac, fiavor 
or color space. This can be done by defining the functional 

eM^^^^ ^j:^;^s^, --' ■ («) 

JVqVqe »^ 



taking the appropriate number of derivatives with respect to external sources 
rji and fjj and putting 77^ = and f]j = [14]. This gives 

n 



where 11 stands for a permutation of the n indices. In particular we obtain for 



the quark condensate (gg) 

m, = {-)iv,c{S,{x,x)}.,=, = (_) ( _;^ ) M / dss^^^ ) (20) 




where X{s) = sA^s)"^ + B^s)"^. fi is the renormahzation scale which we chose 
to be IGeV^. As indicated the trace in eq.(20) is to be taken in Dirac and 
color space, whereas the flavor trace has been separated out. 

Another important consequence from (19) is the fact that the four quark 
operators factorize in the mean field vacuum. For instance in case of T^^^ = 
^ = lfi~2' °^^ finds from (19) and after reapplying the Fierz transformation 
(13) 

(: (g7.f g) (g7.f g) :) = Hi fm' (21) 



i.e. our approach is consistent with the vacuum saturation assumption of 
ref.[l]. The fact, that the four quark condensate factorizes due to eq.(21) is 
independent of the model gluon 2-point function which is used. That means 
that there are no four quark correlations at the mean field level, which is of 
course expected in the large Nc limit. 

Let us now turn to gluonic observables. Because the functional integration 
over the A field in (9) is quadratic for a given quark-quark interaction D we 
can perform the integration over any number of gluon fields analytically. We 
obtain as a generalization of (10) using the same short hand notation 

1.^-1.,.. 1 



jVAAe~¥'''"^+^^= ijD)iy,) eV^^ 
JvAA'e->-'^^^^^ [D{y^,y,)+{jD) 



q\jdj 



(22) 



This means that the gluon vacuum average renders effectively a quark color 
current q'^^-^q together with the gluon 2-point function D. The integration 
over the quark operators is then performed in the mean field vacuum as de- 
scribed above. The remaining expression contains the gluon 2-point function 
D as well as the quark propagator 5*0, which, in turn is given by the self energy 
functions A and B through the Dyson- Schwinger equation (17). The integra- 
tions are most conveniently performed in momentum space. This way we can 



in principle obtain the vacuum expectation value for any gluonic and com- 
bined quark-gluonic operator. It should be stated, however, that the number 
of terms produced by (22) will increase rapidly with the number of gluon fields. 
For instance for the gluon condensate (G^), which contains an A^ integration 
the calculation gets already rather involved. 

In case of the mixed condensate gs{qG^uO'^''q) we have to integrate only over 
powers A^ and A^ and therefore the procedure is still feasible. Applying the 
method described above we find 



-2i)Nc I dz Id'^;^^ glD{z - x)\ii^[SQ{z,x)a^^SQ{x,z)-i^] + 



[+Ai) Nc I dzidz2 glD{zi - x) glD{z2 - x) 

■ ii^[SQ{z2,x)(y^,ySo{x,zi)-i^,So{zi,Z2)'^u] ■ (23) 



In order to calculate this expression we can make explicit use of the Dyson- 
Schwinger equation (15) which determines the mean field vacuum configura- 
tion and which can be cast into the form 



-glD{x-y)S^{x-y) = \\ii,[i:{x - y)] - \\^^tr,[i:{x - y)^,] . (24) 




This eliminates the integration over the gluon 2-point function in favor of the 
quark self energy and therefore strongly simplifies the evaluation of (23). We 
obtain as final result for the mixed condensate in Minkowski space 



9s{qG,.a^^q), = (-) ( ^ ) ^ - y dss- [2A{A - l)s + B'\\ + 

12 1 ds-s" — {2- A)\ . (25) 



In tables 1,2 and 3 we display the result for (gg) and gsiqi^Gq) for three 
different model gluon 2-point functions, which are parametrized by a "width" 
parameter A and a "strength" parameter x [12]. As stated above in all cases 
we fix the pion decay constant in the chiral limit to /,r = 87MeV [13]. The 
renormalization point is in all cases /x = IGeV^. We also show in all three cases 
the values obtained for the chiral low energy coefficients Li, i = 1,3,5,8. As 
we can see the values we obtain for the Li after fixing /jr are compatible with 
the phenomenological ones [13]: 
Li = 0.7 ± 0.5, Ls = -3.6 ± 1.3, L^ = 1.4 ± 0.5, Lg = 0.9 ± 0.3. 



Finally in table 4 we compare with the values which other nonperturbative 
approaches find for (gg) and 5's(gG^i/Cr^'^g): QCD sum rules [3], quenched lat- 
tice QCD [5] and the instanton liquid model [6]. As we can see our results for 
QsiMG^ya^^q) are compatible with the range obtained within these methods, 
whereas the quark condensate {qq) itself is systematically smaller then the 
"standard" value of — (220MeV)^. It should be noted in this context that both 
the QCD sum rules value from ref.[3] and the quenched lattice values from 
ref.[5] refer to /x = IGeV^, whereas the instanton liquid model [6] has a renor- 
malization point which is given by the instanton size and therefore somewhat 
smaller ^ (eOOMeV)^. 



Table 1 

{qq) and QsiqaGq) at ft = IGeV^ for glD{s) 



(47r2d)^, d - 27 



T=. In all cases 



/^ = 87MeV. Also displayed 


are the chira 


low energy coefficients Li,i = 


A 


X 


-im) 


-QsiqaGq) 


Li 


Ls 


L5 


Ls 


[GeV^] 


[GeV] 


[MeV^] 


[MeV^] 


*103 


*103 


*103 


*103 


1 * 10-1 


1.77 


(183)3 


(459)5 


0.79 


-3.76 


2.38 


1.03 


1 * 10^2 


1.33 


(178)3 


(457)5 


0.80 


-4.01 


1.88 


0.93 


1 * 10-4 


0.95 


(175)3 


(456)5 


0.81 


-4.23 


1.48 


0.88 


1 * 10-^ 


0.77 


(171)3 


(452)5 


0.83 


-4.41 


1.19 


0.84 



Table 2 

Tlie same as in tab.l for glD{s) = 37r^ ^ e~A. 



A 


X 


-{m) 


-QsiqaGq) 


Li 


L3 


L5 


^8 


[GeV2] 


[GeV] 


[MeV3] 


[MeV5] 


*103 


*103 


*103 


*103 


0.200 


1.55 


(183)3 


(448)5 


0.81 


-4.07 


1.56 


0.80 


0.020 


1.39 


(167)3 


(431)5 


0.83 


-4.43 


0.96 


0.81 


0.002 


1.23 


(151)3 


(395)5 


0.85 


-4.46 


0.82 


0.93 



To summarize: We have considered a truncation of QCD which gives rise to 
an effective quark-quark interaction whose kernel is given by the model gluon 
2-point function. This gluon 2-point function is chosen to reproduce /,r and it 
also reproduces the physics of the low energy mesonic sector. We have shown 
how to obtain vacuum matrix elements for any quark an gluon operator at 
the mean field level, which is consistent with the large Nq limit if the quark- 
quark interaction is fixed. In particular we have found that any gluon operator 
reduces to an interacting multi-quark operator, whose vacuum expectation 
value can then be calculated in the mean field vacuum. We have applied this 



Table 3 

The same as in tab.l for g^D(s) 



3tt 



2 X 



e A + 



iw^d 



sln[s/A^+e]' 



12 

27 



, A = 200MeV. 



A 


X 


-m 


-gsiqcrGq) 


Li 


Ls 


^5 


^8 


[GeV^] 


[GeV] 


[MeV^] 


[MeV^] 


*103 


*103 


*103 


*103 


0.200 


1.65 


(173)3 


(458)5 


0.81 


-4.03 


1.66 


0.83 


0.020 


1.55 


(160)3 


(448)5 


0.82 


-4.39 


1.13 


0.84 


0.002 


1.45 


(151)3 


(432)5 


0.85 


-4.43 


0.97 


0.88 



Table 4 





-(g-g)^ [MeV] 


-{gsqaGq)'^ [MeV] 


this paper 


150 - 180 


400 - 460 


QCD sum rules [3] 


210 - 230 


375 - 395 


quenched lattice [5] 


225 


402 - 429 


instanton liquid model [6] 


272 


490 



technique to calculate the mixed quark-gluon condensate and found a value 
compatible with the ones of various other nonperturbative approaches. 
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